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Abstract 

In order to analyze data on joint charged-particle/photon distribu- 
tions from an experimental search (T-864, MiniMax) for disoriented 
chiral condensate (DCC) at the Fermilab Tevatron collider, we have 
identified robust observables, ratios of normalized bivariate factorial 
moments, with many desirable properties. These include insensitivity 
to many efficiency corrections and the details of the modeling of the 
primary pion production, and sensitivity to the production of DCC, 
as opposed to the generic, binomial-distribution partition of pions into 
charged and neutral species. The relevant formalism is developed and 



tested in Monte-Carlo simulations of the MiniMax experimental con- 
ditions. 
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1 INTRODUCTION 



There has recently been renewed interest in semiclassical mechanisms of pion 
production in high-energy collisions of hadrons and of heavy ions [1-11]. One 
hypothesis in particular is that pieces of strong-interaction vacuum with an 
unconventional orientation of the chiral order parameter may be produced in 
high energy collisions [12]. This disoriented chiral condensate (DCC) is then 
supposed to decay into a coherent semiclassical pion field having the same 
chiral orientation. 

The primary signature of this mechanism is the presence of large, event- 
by-event fluctuations in the fraction, /, of produced pions that are neu- 
tral. Conventional mechanisms of particle production, including those used 
in standard Monte Carlo simulations, predict that the partition of pions into 
charged and neutral species is governed by a binomial distribution which, in 
the limit of large multiplicity, leads to a sharp value of / ~ 1/3. We refer to 
this as generic pion production. On the other hand, for the decay of a pure 
DCC state the distribution of neutral fraction is very different, following an 
inverse square-root law in the limit of large multiplicity [1-7, 12]. Some other 
production scenarios involving the common feature of coherent final states 
lead to identical / distributions [9, 10, 13-15]. 

Sophisticated phenomenological techniques have been developed in order 
to study the properties of multiparticle final states, and much has been done 
on multiplicity distributions, correlations, and fluctuations [16-20]. Most of 
the practical studies, however, have considered the properties of a single 
species at a time. In the case of DCC, formal tools for the study of the 
joint distribution of neutral and charged pions are required, and here there 
is much less data and corresponding analysis experience [21-25]. 

The authors of this paper comprise the MiniMax collaboration (Fermilab 
T-864), who for the last three years have carried out an exploratory search 
for signals of DCC at the CO area of the Tevatron collider [26]. The heart of 
our detector is a telescope of 24 multiwire proportional chambers (MWPC), 
with a 1.0 radiation-length lead converter inserted after the eighth plane, so 
that charged tracks and converted photons can be counted event-by-event. 
The acceptance in the lego space of pseudorapidity, rj, and azimuthal angle, 
0, is roughly a circle of radius 0.65 centered at r] — 4.1. In 1995-1996, 8 
million triggered events at y/s = 1.8 TeV were recorded. The purpose of this 
paper is not to report the results of this experiment, but rather to describe 
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the techniques we are using as the basis of our data analysis strategy. We 
believe these techniques have much wider applicability and may be of value 
in other searches for DCC signals. 

Even from this very brief description of the experiment, it should be clear 
that we face many challenges in trying to infer either the presence or absence, 
within limits, of DCC signals from the data. These include the following: 

(a) The MiniMax acceptance is small, so that it is improbable that both 
7's from a tt° enter the detector acceptance. 

(b) The conversion efficiency per 7 is about 50%. 

(c) Not all 7's come from 7r°'s. 

(d) Not all charged tracks come from 7r ± 's. 

(e) Because of the small acceptance, the multiplicities are rather low, so 
that statistical fluctuations are very important. 

(f ) Detection efficiencies for charged tracks and 7's are momentum-dependent 
and are not the same. 

(g) Efficiency functions may be dependent upon the observed multiplicity 
or other parameters. 

(h) The efficiency for triggering when no charged track or converted 7 is 
produced within our acceptance is relatively low and different from 
that for events in which at least one charged particle or converted 7 is 
detected. 

Nevertheless, we find that there do exist observables which are robust in 
the sense that, even in the presence of large (uncorrelated) efficiency correc- 
tions and convolutions from produced 7r°'s to observed 7's, the observables 
take very different values for pure DCC and for generic particle production. 
Each such observable is a ratio, collectively referred to as R, of certain bi- 
variate normalized factorial moments, that has many desirable properties, 
including the following: 

1. The i?'s do not depend upon the form of the parent pion multiplicity 
distribution. 
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2. The i?'s are independent of the detection efficiencies for finding charged 
tracks, provided these efficiencies are not correlated with each other or 
with other variables such as total multiplicity or background level. 

3. Some of the i?'s are also independent of the 7 efficiencies in the same 
sense as above. In the remaining cases, the i?'s depend only upon one 
parameter, £, which reflects the relative probability of both photons 
from a tt° being detected in the same event. 

4. In all cases R is independent of the magnitude of the null trigger effi- 
ciency; see comment (h) above. 

5. The ratios R possess definite and very different values for pure generic 
and pure DCC pion production. 

The idealizations implicit in the realization of properties 1-5 include the 
assumptions that particles other than pions can be ignored, that there is no 
misidentification of charged particles with photons, and that the production 
process can be modeled as a two-step process, with a parent-pion multiplicity 
distribution posited, followed by a particular charged/neutral partitioning 
of that population by, e.g., a binomial or DCC distribution function. In 
addition, there is the vital assumption that detection efficiencies for finding 
a 71-^ or 7 do not depend upon the nature of the rest of the event. The validity 
of these idealizations is not contradicted by the simulations presented in this 
paper. This idealized model thus appears to be a good basis for a first-order 
analysis of the properties of the ratios R. We anticipate that this will remain 
true for observations more general than those of the MiniMax experiment. 

The layout of this paper is as follows: In Section 2 we review the conven- 
tional formalism [16-20] of single-variable generating functions and factorial 
moments used in describing global multiplicity distributions. We then de- 
velop the extensions required to describe the bivariate case of distributions 
of 7T ±y s and 7r°'s. The modifications needed to accommodate the decay of 
7T°'s into 7's, as well as the inclusion of less-than-perfect detection efficien- 
cies for charged tracks and 7's, are considered in Section 3. In Section 4 
we introduce the robust observables R and demonstrate their sensitivity to 
charged-particle/photon correlations and their insensitivity to detection in- 
efficiencies and the overall aspects of the primary production process for a 
wide class of production models. The DCC distribution is shown to fall 
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into that class, but with distinctly different values of the i?'s that clearly 
distinguish it from the generic distribution under realistic experimental con- 
ditions. Generalizations of the formalism which allow for the admixture of 
both generic and DCC charged/neutral production are considered in Section 
5. In Section 6 we estimate, by Monte Carlo simulation as well as by use 
of the UA5 charged-particle/photon data at 200 GeV and 900 GeV [25], the 
effects on the i?'s from the realistic complications discussed in the preceding 
paragraph. Concluding remarks are made in Section 7. A number of new re- 
sults concerning the interpretation and representation of the standard DCC 
probability distribution that are needed to establish our results concerning 
DCC production are presented in the Appendix. 



2 GENERATING FUNCTIONS FOR CHARGED- 
PION/NEUTRAL-PION DISTRIBUTIONS 

The entire content of a set of probabilities {P(N)} for the production of 
N particles in a fixed region of phase space can be encapsulated into the 
generating function 

oo 

G(z) = £ z N P(N) (1) 

7V=0 

whose derivatives evaluated at z — 1 yield the factorial moments 

(#G(z)\ 



fi 



dz i / . 

/ z=l 



= (N(N -1)---(N + (2) 



It is often useful to express P{N) as a Poisson transform [27] where one 
introduces a spectral representation in terms of Poisson distributions with a 
weighting function p(p): 

P(N) = I dp p(p) f^e-", (3) 

where 

dp p(p) = 1. (4) 
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The Poisson transform isolates the random statistical fluctuations from the 
physics contained in As an example, the negative binomial parametriza- 
tion 

M = ^Z' 1 (5) 

where X = k/ (N), gives a fairly good two-parameter description of charged 
multiplicity distributions [16, 17]. From (1) and (3) we also obtain a spectral 
representation for the generating function: 

G{z)= / dupMe!*- 1 ), (6) 
Jo 

where now the factor e^ z ~^ reflects the purely random character of the 
Poisson distribution. 

The generating function formalism has been widely used to study charged- 
hadron multiplicity distributions [16-20, 27]. We next generalize this formal- 
ism to bivariate distributions of charged and neutral pions. Among our mo- 
tivations for doing this is the simple manner in which detection inefficiencies 
and particle decays can be handled with generating functions [27]. These fea- 
tures are particularly important in dealing with the MiniMax experimental 
situation. Here the parent 7r°'s are not reconstructed from the observed 7's 
and the efficiencies for detecting both the charged particles and the photons 
are less than perfect. These extensions are taken up in detail in succeeding 
sections. Some earlier work in this connection is contained in Refs. [21-25] 

Let p(n c h, no) denote the probability distribution for the occurrence of n c h 
and no charged and neutral pions, respectively, in a multiparticle event within 
a given phase-space region. As in the single-variable case, the content of 
this bivariate distribution can be conveniently represented by the generating 
function for factorial moments defined by 
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G(z ch ,z ) = V p(n ch ,n )z2 h z%°. (7) 



The partial derivatives of G(z c h, z ) evaluated at z ch — z — 1 generate the 
factorial moments referring to charged (ch) and neutral (0) particles: 

t ( u a\- ( d^Gizc^Zo)} 
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For example, we have, 

/i,o(c/i,0) = (rich), 

/ o ,i(c/i,0) = (n ), 

fi,i(ch,0) = (n ch n ), 

/ 2 ,o(c/i,0) = (n^n^-l)). (9) 

Next, let -P(iV) be the probability for producing a total of N pions with 
any distribution of charge among them. Then p[n ch , n ) can be written as 
the product of two disjoint probability distributions: 

P(n ch , no) = P(N)p(n ch , n ; N), (10) 

where N = n c h + n , and 

oo 

£ P(N) = 1, (ii) 

7V=0 

oo 

E 8N, nch + no p(n ch ,n ;N) = l. (12) 

n c h=0,no=0 

What we call the generic model for the charged-neutral distribution p(n c h, no; N) 
involves no correlations, namely, a binomial (Bin) distribution of n c h and no: 

PBin(n ch ,n ;N) = ( ^ ) T°(l " fT ch - (13) 

Here / is the mean fraction of 7r°'s, which is expected to be about 1/3 as a 
consequence of isospin symmetry. If we substitute (13) into (10) and explic- 
itly denote the dependence on /, the generating function (7) becomes, in the 
binomial case, 

G Bin (z ch ,z ; f) = Y,P(N)ifz + (1 - f)z ch ] N , (14) 

N 

which only depends on the linear combination 

C = fz + (1 - f)z ch . (15) 

Conversely, if a generating function G(z c h, Zo) is a function only of (, the 
charged and neutral pions are binomially distributed. 
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If P(N) is a Poisson distribution, \nGBin(z c h, zo; f) is linear in £. The 
simulations of generic production described in Section 6 yield generating 
functions that, to good approximation, depend only on a fixed linear com- 
bination of z c h and z ; the incorporation of the modeling of the MiniMax 
detector into these simulations is found to alter this linear behavior slightly. 

Much of the simplicity of the generic case is also realized for what can be 
called the binomial transform 

p(n ch ,n ;N) = ( * ^ £ dfp(f)f n °(l - /)"■*, (16) 
of the normalized distribution p(f), 

[ 1 dfp(f) = l. (17) 
Jo 

This leads to a wide class of possible pion factorial-moment generating func- 
tions, namely 

G(z ch ,z )= dfp(f)G Bin (z ch ,z ;f), (18) 
Jo 

where G Bin{z c h-, zq] f) is given by (14) with / replaced by an arbitrary /, 
< / < 1. Combining (3) and (14) we obtain 

G(z ch ,z )= dii pfa) df p(f) e^/Ml, (19) 

JO JO 

where again ((f) is given by (15) with / replaced by an arbitrary /. 

The forms of p(f) and p(ji) depend on the production model and the 
detector. The uncorrelated, generic case (14) corresponds to p(f) = S(f — f), 
where / is some fixed value of /. 

It is shown in the Appendix that for a simple DCC model [1-7] and with 
a sampling prescription appropriate to the experimental situation, p(f) = 
l/(2 v / 7). Although the same bivariate distribution is realized in other hadronic 
production models leading to coherent states [9, 10, 13-15], we refer to this 
case as the DCC model. We note that in the DCC model (no) = 2(n c h), just 
as in the generic case for / = 1/3. 

It is quite possible that the parent pion distribution P(N) or, equivalently, 
p(/x), will be different for the DCC and generic production mechanisms. This 
distinction is important for our considerations of admixtures of the two mech- 
anisms. We investigate some possible scenarios for such admixtures in Section 
5. 
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3 GENERATING FUNCTIONS FOR CHARGED- 
PION/PHOTON DISTRIBUTIONS 



For a detector that is designed to observe charged particles and converted 
7's within its acceptance, events are classified only according to the numbers 
of charged particles and photons, n c h and n 7 , respectively. With sufficiently 
large statistics we can determine probabilities, p(n c h, n 7 ), for observing these 
combinations over some portion or all of the available phase space. 

In order to obtain the charged-pion/photon generating function, incor- 
porating both tt^ 1 and 7 detection efficiencies from G(z ch ,z ), we extend 
Pumplin's cluster theorem [27] to the bivariate case. Consider a generating 
function G(z ch , z ) that refers to charged and neutral "clusters." Suppose, 
for the sake of simplicity, the charged clusters decay in a number of ways into 
charged particles and likewise for the decay of neutral clusters into neutral 
particles. For each of these decay scenarios there is a probability distribution 
and a corresponding generating function, g c h(z c h) or go(zo), respectively. The 
bivariate generating function of the factorial moments of the final charged- 
neutral particle production is then G(g c h(z c h), go(zo)). If the charged clusters 
do not decay, then g c h(z c h) = z c h- On the other hand, 7r° — > 77 with perfect 
photon detection efficiency correponds to go(z 7 ) — z %- 

More realistically, there is a probability e c h for observing a given primary 
charged pion in the detector and a probability 1 — e c h for not observing it. 
These possibilities can be regarded as the two "decay" modes of the primary 
charged pion which is otherwise regarded as stable. Similarly, there are 
probabilities e m , m = 0, 1, 2, with 



for observing m photons from a 7r° decay and each possibility can be re- 
garded a decay mode of the 7r° cluster. If these probabilities are identified 
with what we suppose are the independent, i.e., uncorrelated, efficiencies for 
the respective detection options, the generating function for the distribution 
of observed particles, including efficiencies, is obtained from G(z c h, z ) by 
replacing z c h by the generating function 



e + ei + e 2 = 1, 



(20) 



g C h{z c h) — (1 — Ccfe) + tchZch, 



(21) 
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and zq by the generating function 

goizj) = e + eiz 7 + e 2 z\ (22) 

For the class of production models characterized by (18), the preceding 
considerations lead to the following factorial-moment generating function for 
the distribution of observed charged pions and photons: 

G obs (z ch ,z 7 ) = ( dfp(f)G Bi n(gch(z ch ),g (z y ); f). (23) 

J 

The charged-pion/photon factorial moments are 

d^G{z ch)Zl ) \ 

which introduces the bivariate indexing with respect to charged particles 
and photons employed henceforth. For example, the two lowest orders of 
factorial moments are: 



f lfl (ch n ) = (nch) = (1 - f)e ch (N), (25) 

/ 0> i(cfc,7) = W = (/)(e 1 + 2e 2 )(7V), (26) 

f 2 ,o(ch n ) = (n ch (n ch - 1)) = ((1 - f) 2 )el h (N(N - 1)), (27) 

h,\{ch,l) = (non^) = (/(l - /))e ch (e 1 + 2e 2 )(N(N - 1)), (28) 

/o, 2 (c/i, 7 ) = KK - 1)) = (/ 2 )(ei + 2e 2 ) 2 (N(N - 1)) + 2e 2 (/)(7V). (29) 

In Eqs. (25)- (29) the overall statistical averages for the charged, photon, and 
charged-photon factorial moments are expressed, in an obvious notation, in 
terms of the independent moments taken with respect to the P(N) and p(f) 
distributions. 

Finally, we turn to the effect of the MiniMax trigger on these considera- 
tions. The MiniMax trigger requires, among other things, a coincidence in the 
signals from scintillator counters located behind both the converter and the 
entire tracking telescope. In consequence, events in which no charged particle 
or converted 7 goes through the acceptance of the detector are triggered with 
different (and lower) efficiency, e, than events in which either a charged par- 
ticle or 7 conversion products go through the aperture. An effective model 
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for the effect of the MiniMax trigger on the probability, p obs (n c h, n 7 ), for 
observing an event with n c h charged particles and n 7 converted 7's passing 
through the acceptance is given by the proportionalities 

p*"9(0, 0) = eap obs (0, 0), n ch = n 7 = 0, (30) 

and 

p tri9 (n ch , n y ) = ap obs (n ch , n 7 ), n ch + n 7 > 0. (31) 

Here p tri9 (n c h,n y ) is the measured probability of seeing an event, including 
the effects of both the trigger and the particle detection efficiencies, while 
p obs (n c h, n 7 ) presumes perfect triggering. If 

a=[l + (l-e)/ te (0 ) 0)]- 1 ) (32) 

ptng w jjj De p r0 p er ly normalized if p obs is. 

The bivariate factorial moments transform homogeneously under the trans- 
formation (30)- (32) incorporating differential trigger efficiencies, 

fi,j(ch, 7) -»• afij(ch, 7)- (33) 



4 ROBUST OBSERVABLES 

The second-order factorial moments (25)- (29) represent the lowest-order cor- 
relative effects among charged pions and photons. We see from (29) that 
the gamma-gamma correlations are distinguished by the term 262(f) (N) for 
observing the two photons from a single neutral pion, so that this average 
will not be a component of a robust measure involving only first and second 
order moments. This suggests the construction of a measure from the mo- 
ments (25)- (28) in the form of a ratio in order to cancel out as many effects 
as possible, apart from the p(f) averages, that reflect the particular details 
of the production mechanism. 
Consider, then, the ratio 

n 1 = ( n chn~/)(n ch ) ^ 
(n c h(n c h ~ !))(%) ' 
For generating functions of the form (23), we find from (25)-(28) that 

r (/(I -/)>((!-/)> m) 

Xl ~ ((i-m(f) ' (35) 
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an expression in which all reference to the background distribution P{N) and 
the efficiencies ei, 62, and e c h have cancelled out. Further, we see that 

n,i -> n,i (36) 

under the transformation (30)-(32) so that r^i is a "robust observable" in 
the sense referred to in Sec. 1. 
It follows from (35) that 

ri,i < 1, (37) 
where the equality is realized for generic pion production, p(f) = S(f — f), 

r^i(generic) = 1, (38) 

independently of /. The realization of the limit (38) in the UA5 data at 
200 GeV and 900 GeV [25], and in Monte Carlo simulations at 1.8 TeV, 
both of which include nonpionic sources of charged particles and photons, is 
considered in Sec. 6. 

For a DCC distribution, p(f) = l/(2v7), one finds 

r hl (DCC) = \. (39) 

This clearly distinguishes the pure DCC and generic distributions. 

The values (38) and (39) represent the limiting extremes of a mixture 
of generic and DCC distributions. Generally, broad (DCC) and narrow 
(generic) statistical distributions can be distinguished in a mixture of the 
two by means of higher-order moments that are sensitive to the tail of the 
charged-particle/photon distribution. Robust combinations of these higher- 
order moments that are generalizations of r\ i will be of greatest practical 
value in an analysis of data in which a discernable fraction of DCC form is 
expected to appear. 

Let us first note that the normalized factorial moments 

(Af-i+l)) 

(Ny v ; 

are unity if the parent distribution P{N) is Poisson. Therefore, deviations 
from purely random fluctuations are measured by the departure of the Fj's 
from unity. A bivariate generalization of the F^s is given by 

_ (n ch (n ch - 1) . . . (n ch - i + 1) n 7 (n 7 - 1) . . . (n 7 - j + 1)) 

{rich} (n^y 
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In particular, one finds that 



F,o - 1(rrjr (42) 



^ = * y,: (43) 



and 

jj+i - fY ) 

(/)((! -/)>* 

where Fj refers to the ith normalized factorial moment (40) of the P(N) 
distribution for the total multiplicity. We note that 



/••„ - ; '/•;, (44) 

under the transformation (30)- (32). 

Evidently, r± t i = Fi ; i/F 2i0 . From (42) and (43) we find a generalization 
of ri,! to a family, R, of robust observables: 

M (/) ((1 " fY +1 ) ' 1 J 

Moreover, one finds that for all i > 1 

r iA (generic) = 1, 
nADCC) = (46) 

in the two cases. Thus, r^i becomes more sensitive to the difference between 
DCC and generic production mechanisms with increasing order of the mo- 
ments. This reflects the broadness characteristic of the DCC distribution in 
the neutral fraction / compared to the generic case. 
The ratios 

rij = tP^ (47) 

■ri+3,0 

are not robust because the moments Rij for arbitrary % and j are not inde- 
pendent of the photon detection efficiencies. However, the terms involving 
these efficiencies can be expessed in terms of only one combination of these 
parameters, namely 

^ = ( £l + 2e 2 ) <n 7 >' (48) 
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along with the mean number of photons, as 



m ((1 - fYp-" >) 

^o' J '"^ ' ' "'((1 -/))''</)•• 



The coefficients obtained from the identity, true for any differentiate 

function, D(z 2 ), 

dW(z 2 ) _ om /^ y -2 m ^ m D(z 2 ) 



The first few Cj jTn are 

c i,o = 1, 

c hl = j(j -l)/2, 

c j:2 = 3j!/4!(j-4)!. (51) 

One can use the ratios r it /s in the analysis of experimental distributions, 
with the understanding that the parameter £ is to be determined from the 
data. Generally, we have the bounds and limiting values 

r\ j {generic) > 1, (52) 

(r id (generic) )^ =0 = 1, (53) 

and 

(r -(iX7C7))^ = ,!( ^" 1)!! . (54) 

5 SENSITIVITY TO DCC ADMIXTURES 

We next turn to the question of what can be said about robust observables 
when there is an admixture of DCC and generic multipion production. There 
is considerable theoretical uncertainty about how such an admixture would 
arise in hadronic collisions and so there are many possibilities for extending 
the development given in the preceding sections. Our objective in this sec- 
tion is only to provide a formalism in which the sensitivity of experimental 
results to the presence of DCC or some other anomalous mechanism can 
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be investigated. Thus, it will suffice to address this question only in the 
context of a few simple limiting models of pion production containing both 
generic and DCC components. Specifically, we consider modifications of the 
generating-function formalism we have developed in the preceding sections 
in three different scenarios for mixing DCC and generic multiparticle pro- 
duction. Then we examine the impact of these modifications on the values 
of the robust observables. 



5.1 Exclusive Production 

First, let us consider the possibility of what we refer to as exclusive produc- 
tion. That is, in a given event, particle production is either the result of 
the formation of a DCC with probability A, or it is generic with binomially 
distributed charged and neutral particles with probability 1 — A. The picture 
of exclusive production could be regarded as a first-order phenomenology of 
very high-energy cosmic-ray interactions, which seem to divide themselves 
into what appear to be generic and anomalous classes [29]. 

The generating function for the exclusive production of charged pions 
and the photons resulting from it decay is simply the weighted sum of the 
generic and DCC generating functions: 

Gexcl( z ch, z 7 , A) = (1 — X)Gg eneric {Z c h, z y ) + \Gocc( z ch, z y ). (55) 
Here G generic (z c h, z 7 ) and Gocc( z ch, z 7 ) are obtained from (23) for the cases 
p{f) = $(f ~~ f) an d p(f) = V(2\/7)> respectively, and where the distribu- 
tions P{N) of the total number of pions are generally different in the two 
cases. 

The expressions for the moments r^i obtained using G exc i(z c h, z 7 , A) in- 
terpolate between the generic and DCC limits as A varies between and 1. 
For example, since 

m d = n \).ir:' ■ \n:r r 

fDCC 

= CT + ^-l)), (56) 
it follows, using the results of Section 3, that one can write 



r excl ( ^ + A[ 15/(1-/) (N(N-l))^ + A l 3 (l-/) L Ji ( , 7) 

1,1 n i \( 8 mN-i)) Dcc iNiri | u i m DCC m 1 ' 
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Note that this expression explicitly depends on the relative size of the DCC 
and the generic factorial moments. Technically, this ratio is no longer "ro- 
bust" in the sense of the preceding section. However, it still does not depend 
upon efficiency corrections. In addition, the extra dependence will be an 
advantage if DCC dominates the high-multiplicity tail of the distribution. 



5.2 Independent Production 

A second possible production scenario is where the occurrence of DCC in an 
event is independent of the pions that are produced generically. Independent 
production implies that the probability Pdcc(N) for producing N DCC pions 
is independent of the probability P ge neric(N) for producing iV binomially 
distributed pions, so that the generating function factors into a product, 

(zch,Zy)G D cc(zch,Zy). (58) 

Thus, we find 

fit = i t ( 1 ) ( I ) /S,-X? c m 

Hence, using the results of the previous sections, it follows that, for example, 

{N) Gen {N) DCC 2 ± 2{N(N-1)) DCC 1H 2{N) DCC 1 

ind = [ (N(N-l))^ ^3(1-/V 3f' _ 15/(l-/)(iV(iV-l))Q e " J[ 3(l-/)(jV)g«" J 

1,1 r (tyGen^DCC , 4 \ | 8{N(N-1)) DCC 1M 1(N) DCC 1 

(60) 

Again the sensitivity to the independent production of DCC is dependent on 
the ratios of DCC and generic factorial moments, but not to the efficiency 
corrections. 

We note that in the independent production model 

(z c h, Zj) + lnG DC c(zch, Zy), (61) 

which suggests an analysis in terms of a bivariate generalization of single- 
variable cumulant moments [16, 18-20]. We define bivariate cumulants for 
i + j > by 

•*-(& hG L^ (62) 
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/From (58) we see that in this production scenario, the cumulants are addi- 
tive: 

k™f = k°™ eric + k?f c . (63) 

For single-variable probability distributions, cumulants reflect non-random 
correlations in that they vanish for a Poisson distribution. In the bivariate 
case their properties as a measure of correlations are not so direct. 

As with the bivariate normalized factorial moments (41), we introduce 
normalized bivariate cumulant moments: 

Kij = (n ch r (n,r [jfy lnG ) • ( 64 ) 



In the independent model we obtain for Kf^ the weighted sum 

Ktf = K h ^K CC + (1 - WU - \) j Kr TiC i (65) 

where 

are the fractions of the mean charged or photon multiplicities attributed to 
the DCC. 

The formulae for the normalized cumulant moments for DCC and generic 
subsamples are obtained in a straightforward manner. As before, most of 
the efficiency corrections cancel out. However, the cumulant moments do not 
scale homogeneously under the differential trigger inefficiency characteristic 
of MiniMax. While this is disadvantageous for the early MiniMax analyses, 
there is reason to expect that they will be eventually of substantial utility in 
MiniMax as well as in other experiments. 



5.3 Associated Production 

A third possiblity for the contamination of a DCC signal by generic multipar- 
ticle production is what can be called associated production. For example, 
in the Baked Alaska model [8] the number of DCC pions is estimated to scale 

as 

N DCC ~ (AW-) 3/2 - (67) 
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A simpler case, which is also a credible scenario, is where the amount of 
DCC production is, on the average, proportional to the amount of generic 
production. It then follows using the cluster theorem [27], that 

G assoc (z ch ,z 7 ;X) = f df b p b (f b ) f df d p d (f d ) 

Jo Jo 

oo 

]T P(N)[(1 - X)g b (z ch ,z 7 )+Xg d (z ch ,z 7 )] N , (68) 



N=0 



where, 



g A (z c h, Zj) = fAgo(z-y) + (1 - fA)g c h(z ch ), (69) 

p b (f b ) = 5(f b -f), (70) 

Pd (f d ) = 1/(2^), (71) 

and the index A takes the values b and d in the binomial and DCC cases, 
respectively. As before, one can carry out the calculation of the robust ob- 
servables which results in formulae that interpolate as a function of the frac- 
tion, A, of DCC admixture between the generic and DCC limits. Note that 
in this case there would be only a single parent P(N), common to both the 
generic and DCC production. Using the results of the previous section, one 
can calculate, for example, 

rasS o C(x) = [(1 - A) 2 /(l - f) + |A(1 - A)(l + f) + ^A 2 ][(l - A)(l - f) + IX] 
[(1 - A) 2 (l - /> + |A(1 - A)(l - /) + Ap][(i _ \)f + | A ] 

(72) 

which, in contrast to the other two cases, (57) and (60), is a fully robust 
observable. 



5.4 Other Particles 

A similar framework can be used to discuss the sensitivity of the predictions 
to the production of particles other than pions. This is of potential concern, 
since K and i]° production may be a substantial fraction of pion production 
[17, 25, 30]. In particular, the r]°/ii ratio can be quite large leading to an 
excess of gammas over the case of pions alone, where (n 7 ) = (n c h)- 

Relatively little is known about K and i]° distributions at the highest 
energies, especially in forward directions, so, while an independent produc- 
tion model might be more accurate, we will limit our considerations at the 
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moment to the context of an "associated" production model. In essence, we 
are thus assuming that a system of parent partons is created in the collision 
process, and that this system then evolves into a system of N hadrons with 
probability P(N), with the hadrons independently partitioned into various 
species. 

Let the index % run over the various types of hadrons that are produced. 
The i'th type of hadron is produced with relative probability A« (with J2i K — 
1). These hadrons then decay into charged particles and 7's, and each species 
of hadron is characterized by a generating function for detecting the products 
of that species, 

9i (z ch , zj = £ £ eg, >B7 *3* z^, (73) 

where <?i(l, 1) = 1. Then the observed generating function, neglecting DCC 
production, can be written as 

G obs (zch,z 7 ) = £P(7V)[£A^(^,z 7 )] Ar . (74) 

N i 

We can now make a few observations about the impact of contamination 
of the predictions that arise from K and 77 production. The following esti- 
mates of the effects of various particle types on the magnitude of r^i draw 
upon the simulations specific to MiniMax reported in Section 6. 

First, we note that the K ±, s, which are seen simply as charged particles in 
MiniMax, appear just as another source of charged particles from the collision 
point and so modify the neutral fraction, but are otherwise benign. Similarly, 
the Kls have a decay length much longer on average than the length of the 
MiniMax detector. In consequence they are only detected, but not identified, 
when they interact strongly in the converter used to identify photons. On 
the relatively rare occasions when K^s do interact in the converter, they 
are misidentified as 7's. This will also influence the net neutral fraction 
that is observed, but is also otherwise benign. In conclusion, the associated 
production of K ±, s and K^s will not change the values of the r^j's predicted 
in Section 4 for "generic" production. 

The case of K$ production is rather interesting since the K$ decay modes, 
K s — > 7r + 7r~(69%), K s — > 7r°7r°(31%), are essentially those of an isosinglet 
DCC with one pair of pions. That is, in regard to the statistics of the 
particles produced, K s decays are essentially identical to those of the smallest 
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conceivable domain of DCC's. As such, Ks production is in principle of 
interest from the point of sensitivity to very small domains of DCC. 

Let us consider associated production of K^s with fraction A^ s . The 
generating function for studying the modification of generic production is 
thus 

Gr s ( z c, z^; \ Ks ) =^P(N)[(l — \ Ks )g gen (z c , z^) + \ Ks g Ks (z c , z y )] N . (75) 

N 

Using previous methods, one finds that 

(r^fX))- 1 = 1 + 2(N)\ Ks e%$ 

(N(N - 1)>((1 - X Ks ){\ - f)e ch + X Ks (e^ + 2eg))= 

(76) 

which is manifestly not robust. 

K$'s are not DCC domains; they are, rather, particles of well-defined 
mass and a lifetime such that most of them have decayed before reaching the 
MiniMax detector, and their decay products have strong correlations and 
are not vertexed to the collison point. As a consequence, in MiniMax the 
acceptance for 2 charged pions from a single K s is about 4%. Consequently, 
the impact of K s production on the MiniMax systematics is expected to be 
quite small. 

One can similarly study the impact of rf production on the idealized 
predictions of Section 4. The rf has a wider variety of decay modes and 
all of the charged particles and 7's from the decays are collision vertexed. 
Thus g v o(z c , z 7 ) is more complicated, but the calculations follow closely those 
outlined for Ks decays. In addition to having decay modes with more than 
a single charged particle, there are decay modes with intrinsic charged-7 
correlations, as well as the charged-charged correlations which entered into 
the Ks analysis. The conclusion is, nevertheless, much the same. 



5.5 Detector Effects 

Finally, we note that the formalism we have developed can be extended to 
consider contamination due to detector-related effects. For example, in de- 
tectors which identify gamma rays by electromagnetic calorimetry, charged 
hadrons can also be identified as photons when they interact strongly in the 
calorimeter. For example, in WA98 [31], a heavy- ion experiment at CERN 
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which has instituted a DCC search, this is expected to occur approximately 
20% of the time. Such misidentifications can be handled by using an appro- 
priate form of the generating function gi(z c h, z y ). For example, 

fi , 7r ± (^cfe ) Zy) = e 00 + £\fiZ c h + e 1 ^z c f l z~ / (77) 

would be suitable if some fraction eff of the charged pions were tagged as 
both charged particles and photons because of the calorimeter's response. 



6 ROBUST OBSERVABLES IN PRACTICE 

We now turn to the utilization of the robust observables for analyzing collider 
data, both actual or simulated. As we saw in the last section, the assumptions 
made earlier are idealizations that are violated by some types of production 
mechanisms and by less than ideal detector performance. In this section we 
examine the properties of the robust moments in the context of the UA5 data 
and Monte-Carlo simulations of the MiniMax detector in order to assess the 
importance of these violations in practice. 



6.1 7*1,1 from UA5 

For collider energies of 200 GeV and 900 GeV, UA5 measured the inclu- 
sive charged-particle and photon dN/dr] distributions, as well as the cor- 
responding charged- charged and the charged-photon correlation functions, 
C ch , ch (r} ch = 0,r] ch ) and C ltCh (r] 7 = 0,r) ch ), respectively [25]. Here, r) ch and ?? 7 
denote the charged-particle and photon pseudorapidities, respectively. The 
measurements were carried out over about four units of \r) c h\- The mean 
values {n c h) and (n 7 ) can be calculated for different pseudorapidity bins 
using the experimental dN/dr] distributions. Under the assumption that 
Cch,ch(Vii V2) and C 7;C ^(?7i, 772) depend only on the absolute value of \rji — 772 j , 
the second-order moments that enter into r^i can be also be calculated for 
corresponding pseudorapidity bins. Despite large uncertainties in the UA5 
photon data and the validity of our assumptions about the correlation func- 
tions, we find r^i = 1.0 ± 0.10 for the different energies and various bin 
choices. 
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6.2 Simulations 



While we believe the robust observables will find general application in exper- 
imental searches for DCC, we are motivated here primarily by the MiniMax 
experimental situation. In this context, in order to make a rough check of 
the validity of the assumptions we have made in the opening sections, we 
next describe a series of complete simulations of the MiniMax experiment. 

Minimum bias events are generated in PYTHIA version 5.702 and JET- 
SET 7.401 [32]. The output of PYTHIA is then used as input to the simula- 
tion of the detector response using GEANT, version 3.21 [33]. The GEANT 
output is then put through a full tracking and analysis chain. The resulting 
frequency distributions for observing n c h charged tracks and n 7 converted 
photons are then used to calculate the various robust observables. Similar 
studies, in which the output of PYTHIA is replaced or augmented by the 
output of a DCC generator, are also carried out. We find the results of these 
simulations to be in agreement with expectations from our calculations in 
the previous sections. 

6.2.1 Standard Monte Carlo 

PYTHIA is used to simulate the minimum bias collisions at y/s =1.8 TeV. 
Default values are taken for all parameters except that particles with a mean 
decay length greater than 1 cm were not allowed to decay. 

There are no published data on multiparticle production at 1.8 TeV in 
the pseudorapidity interval covered by MiniMax, so there is no independent 
check on the accuracy of the simulations. For recent measurements at 630 
GeV [34], the agreement between PYTHIA and the dN/dr] data, in a range 
of pseudorapidity including that of MiniMax, is less than ideal. Nonetheless, 
the PYTHIA output represents a useful benchmark. 

The particles generated in a simulated collision are then taken as input 
into a GEANT simulation of the detector and its environment. The exper- 
imental data give evidence of a large background of particles aising from 
interactions in material immediately surrounding the detector. Therefore, 
many nearby objects are included in the simulation. GEANT propagates the 
particles through the detector and its surroundings and produces a simulation 
of the data that are produced by the actual detector. Despite care in includ- 
ing all relevant aspects of the detector and its environment, the GEANT 
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data show a smaller number of reporting wires in the MWPC's than do the 
actual data by a factor of two. 

GEANT data are written to a file that is used as input to the same code 
that is used for the analysis of the actual MiniMax data. The analysis pro- 
ceeds in two stages. First, a tracking code is used to find track segments in 
front of (heads) and behind (tails) the converter plane. The output of this 
calculation is used by a second code (vertexer) that determines the number 
of charged particles and 7's observed in the event. In so doing, it counts 
a charged track to be a head that can be joined to at least one tail. A 7 
conversion is taken to be one or more tails emanating from the same point 
in the converter without an accompanying head. Candidate charged and 7 
tracks are required to point to within some given distance from the collision 
point in order to remove secondary particles from material adjacent to the 
detector and fake tracks arising from chance combinations of random report- 
ing wires. The parameters used in the vertexer are determined by optimizing 
the reconstruction of the events generated by PYTHIA and GEANT. 

This track- reconstruction procedure is still under development. It does 
not satisfy all of the assumptions made in Section 1 regarding tracking effi- 
ciency In particular, the reconstruction efficiency may depend on the mul- 
tiplicity and proximity of tracks. 

6.2.2 DCC generator 

DCC production is modelled according to the l/(2y/f) distribution. For the 
present simulation, the DCC domain size in rj — <fi space is taken to be on 
the order of the detector acceptance. The cm. momentum of the DCC is 
directed at the center of the acceptance with a reasonably large p T . We 
assume that the number of pions in the DCC is independent of the central 
pseudorapidity of the DCC. The ratio of the mean energy density of DCC 
pions to that of generically produced pions is then approximately constant; 
we take the ratio to be unity. 

DCC's are generated using what could be called a "snowball" model in 
reference to the low pion momenta in the DCC cm. The number, N^cci of 
DCC pions is chosen using a Poission distribution with mean fiocc- 

The neutral fraction is generated using the transformation method, where, 
if a; is a uniform deviate, / = x 2 is distributed according to l/(2y/J). A 
uniform deviate Xj is then generated for each of the N DC c pions; if Xj < /, 
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the pion is defined to be neutral, otherwise it is defined to be charged. This 
procedure implements the l/(2y/J) distribution exactly; if one takes the 
viewpoint that the isosinglet distribution is more fundamental, then this 
procedure can be viewed as an approximation to it which is valid in the limit 
that the total number of pions is large, and one is sampling a subset of the 
DCC. The actual distribution is, of course, an experimental question. 

Each of the pions is assigned a 3-momentum in the DCC cm. system by 
drawing from a zero- mean Gaussian distribution with a variance (p-p) = Sa^. 

The DCC is then boosted such that the momentum of the DCC cm. is 
in the direction of the center of the MiniMax detector at rj = 4.1, and so 
that the DCC pions have (p?) ~ cr p . If the pions are not too relativistic in 
the DCC cm. frame, the boosted DCC domain is approximately circular in 
i] — (f> space, with radius Rdcc ~ & p Ipt- 

The results we report next are based on Monte-Carlo simulations in which 
<7 p = 0.1 GeV and = 0.14 GeV; hence Rdcc ~ 0.7, the typical radius of 
a hadronic jet. The Poisson mean for the number of DCC pions has the 
value fiDcc = 5.5, which corresponds to an energy density in lego space 
comparable to that of generic production. The Monte-Carlo simulation of 
DCC production is used to generate pure DCC events. These events are then 
run through the same GEANT simulation as the PYTHIA events, except that 
the trigger is not used since no particles go in the p direction. 

6.3 Results 

Once the number of charged tracks and 7's passing into the acceptance is 
determined, the moments and are calculated. Statistical errors are esti- 
mated assuming Poisson fluctuations and the standard propagation of errors 
formalism [35]. 

The results obtained for approximately 5 x 10 4 PYTHIA events and 2 x 10 4 
pure DCC events are shown in Table 1. For purposes of comparison, the 
predicted values for idealized binomial and DCC distributions are included. 
For those ratios involving higher-order moments of the number of observed, 
converted 7's, the predictions are nonrobust, as discussed in Section 4, and 
depend on £, which is determined from the relationship between f 2 , f2,o, 
and <n 7 >, assuming a binomial distribution. The same values are used in 
correcting the DCC predictions for the higher order moments. In particular, 
it is assumed that 2e 2 /(e 1 + 2e 2 ) « 0.08 ± 0.01 obtained from PYTHIA for 
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generic production has the same value for DCC production. This is certainly 
violated in practice, for the simulated DCC pions have significantly lower 
(jp T ) ~ cr p than those generated by PYTHIA, and hence the probablility of 
both 7's from a 7r° decay being in the acceptance, which is reflected in 62, will 
be different. In addition, the F^s are also taken to be the same in the DCC 
case as in the PYTHIA case, which is also clearly a poor assumption. We 
have chosen to display the data in the manner shown, however, in order to 
illustrate the problems which will arise in DCC searches using these moments. 

There is general agreement between the "predictions" based on the analy- 
sis in Section 4 of this paper, and the results of these full simulations. One of 
the striking features of these results is how well the PYTHIA/GEANT sim- 
ulation, which includes resonance production, simulations of detector effects, 
among other features, matches the predictions of a simple binomial model. 

In order to illustrate the effect of an admixture of DCC with generic 
events, where the amount of DCC produced is independent of the amount 
of generic production, DCC domains from the DCC-generator/GEANT are 
added to various fractions of random PYTHIA/GEANT events. This rep- 
resents a mixture of the independent and exclusive models considered in 
Section 5. The effect on the r^i is shown in Table 2. 

These simulations support the expectation that the robust observables 
introduced in this paper will be a useful analysis tool, even though all of the 
technical requirements of robustness may not be met. Thus, these observables 
provide a well-defined framework for describing correlations in such way that 
many systematic uncertainties cancel out. 

7 CONCLUDING REMARKS 

Most of the experimental analyses and theoretical studies of multihadron 
production have concentrated only on charged-hadron production, for which 
the bulk of the data have been taken; for exceptions to this, see [21-25]. The 
questions we have addressed concerning the neutral-hadron component of 
multiparticle production have received little attention, but are vital for our 
MiniMax experiment. 

The robust observables R which are here proposed appear, on the basis 
of the analytic calculations and Monte Carlo simulations we have presented, 
to be of considerable value in all future analyses of charged-particle/photon 
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distributions in high-energy hadron and heavy-ion collisions, and especially 
with respect to the search for disoriented chiral condensate. 

While these observables are manifestly robust, there are still clear limita- 
tions to their use which must be eventually be addressed. We have said little 
about momentum-dependent efficiencies; this will, at the formal level, require 
generating functions to be generalized to generating functionals [22-24]. At 
this level, even the choice of the parent generating functional may have con- 
siderable ambiguity due to a lack of consensus on the underlying physics; 
e.g., can the Poisson-transform structure of Eq. (6) be simply generalized? 

At a more practical level, the issue of correlated efficiencies, especially 
with respect to total multiplicity and background level, is vital. Here the 
features of the individual experiment and its environment are essential, and a 
strong interplay between simulations and the analysis of real data is required. 

Finally, in experiments with large acceptance, even for pure DCC pro- 
duction the chiral order parameter may be different in different portions of 
the r) — <p, or lego, phase space. In this case the formalism we have presented 
must undergo further generalization. 

Nevertheless, we believe that the analysis strategy we have described can 
serve as a very useful starting point for the experimental search for disori- 
ented chiral condensates. 
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APPENDIX: DCC DISTRIBUTIONS 



The distribution 
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where N and n < N are nonnegative integers, was discovered by Horn and 
Silver [15] in the context of coherent-state production models. For this reason 
we will refer to it as the coherent distribution. It was later found that the 
coherent distribution was an appropriate final state for a simple model of a 
zero isospin DCC [6]. In both physical contexts, the distribution is relevant to 
the case of an even total number, 2N, of pions and, necesarily, because of zero 
isospin, to an even number, 2n, of 7r°'s. In the mathematical considerations 
that follow, n and iV are regarded as arbitrary nonnegative integers. 
In [6] it was shown that 

as N, n — > oo, with n/N = f constant, in agreement with the classical 
expectations for a DCC [1-5, 7, 11]. Generally, a bivariate distribution can 
be expressed as a continuous binomial distribution weighted over the infinite- 
sampling limit: 

P[n- N] = ( N n ) jf 1 ^Lf\l - /)<"-»>. (A3) 

The representation (A3) can be established by rewriting (Al) in terms of the 
beta function, 

as 

P[n; N] = ±[* y jB(n + ±l + N-n), (A5) 
which is found using the identity 

The standard integral representation of the beta function [36] yields (A3). 
The identity (A5) establishes the connection between the coherent-state pro- 
duction model of Martinis, et al., [9], for 1 = and the analysis of [6] and 
[15]. The continuous binomial distribution (A3) allows one to calculate all 
averages in the same explicit manner as for the binomial distribution for a 
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particular / and then integrate the result over / with the indicated weighting 
leading to exact results for the various moments. The direct use of (Al) to 
calculate averages is awkward. 

A problem arises with the original interpretation of distribution (Al) in 
connection with a realistic dectector, or, equivalently, a sampling consisting 
of a finite number of pions. The limited sampling of such a detector means 
that typically one sees only a portion of the particular group of the correlated 
pions that are thought to be the earmark of a DCC. Within that sampling we 
need to find the distribution induced by the DCC and with it we can carry 
out a generating-function analysis. We show that the coherent distribution 
is self-similar in that the neutral/charged distribution of a finite number of 
pions chosen from a sampling space distributed using the limit of (Al) for 
N — > oo, is given in fact by (Al), but now with iV and n regarded as the 
total number and the number of 7r°'s, respectively, whether or not they are 
even or odd. 

In support of these remarks, let us consider the problem of the neu- 
tral/charged distribution of an arbitrary subset, even or odd, of a DCC cor- 
responding to 2N pions that are distributed according to P[n; N]. Suppose 
then that because of limited sampling we observe n t < 2N pions. The joint 
probability distribution function for finding no neutral pions and n c h = fit— no 
charged pion's is then a product with P[n, N] of the hypergeometric distri- 
bution [35] of the two relevant binomial samplings: 

N-\n c 

Q[n ;n t ;N} = ^ 

n>in 

where realizing equality in either of the limits is possible only when these 
limits are even. The nature of the summation limits in (A7) complicates a 
direct proof of the correct normalization, viz., 

nt 

J2 Q[n ;n t ;N} = 1, (A8) 

no=0 

however, (A8) has been verified numerically. 

Because of the limited sample one cannot regard N in (A7) as known. 
Therefore the case where all that is known is that N 3> 1 is of special interest. 





P[n;N], (A7) 
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In this case we find using (A2), the Stirling approximation, and passing to 
the continuum limit of /, that 



Q[n ;n t ;N] -> P[n ;n t ], 



(A9) 



where 



P[n ; n t ] 



(2)I»$7 ril - f) "*~"° K {AW) 



which has precisely the same form as (A3). Here, however, the respective 
functional parameters are the number of neutral and total pions sampled from 
the DCC, rather than half those numbers as they are for all of the pions 
of a full DCC. Thus the induced representation (A10) is a quasi-coherent 
distribution that goes over to the classical DCC distribution (A2) in in the 
infinite- n t limit, which in practice may not be too large, due to the accuracy 
of the Stirling approximation for fairly small numbers. 

The similar forms of Eqs. (A3) and (A10) shows that, in regard to an infi- 
nite sampling space, the coherent distribution generates a self-similar induced 
distribution. In addition, the procedure used to arrive at (A10) indicates how 
one uses the continuum limit of the coherent distribution to define a sampling 
of a finite number of pions from an infinite sampling space. This remark then 
also explains the use of the form (A3) when it is applied to the full DCC: 
It represents a sampling algorithm carried out by means of neutral pairs of 
pions to induce a DCC of a finite, even number of pions out of the infinite 
sample. 

The distribution (A10) refers to a collection of pions that need not have 
net zero charge, the signal characteristic of a full DCC, but yet makes no 
reference to the total charge. For the sampling algorithm used to obtain 
(A10), the absolute magnitude of total charge will obviously be binomially 
distributed about zero if n c h is even, and unity if it is not; this extended form 
of (A10) should be used when the sign of the pions can be distinguished. 
When they cannot, the means and variances have interpretations that are 
different from a DCC. 

Finally, let us weight P[no] n t ] with respect to a parent distribution P[n t }. 
Then the relevant generating function is 



GDCc(Zch, Zq) 



E K,n ch+ n Q P[n ; n t }P[n t }z2 h z%°. 



(All) 



n t ,n ch ,n =0 
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Representation (A10) when combined with (A7) yields 



GdCc(Zc!i, Zq) — / r-p GBidyZchi Zo] /); (^12) 

Jo 2V/ 

which we interpret as the generating function of the factorial moments of 
the numbers of charged and neutral pions sampled from a very large DCC 
sample space. 

The distribution p(f) = l/(2y/J) has been associated with the decay of a 
DCC in the classical limit. Thus, the generating function (A12) can be con- 
sidered applicable to the situation in which the phase-space domain of the 
particles resulting from the DCC is very much larger than the acceptance 
of the detector. Then one can picture DCC production as corresponding to 
an event distribution for which the neutral fraction / is a random variable 
distributed according to 1/(2a/J), a depiction reflected in (A12). 
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Table 1. Robust observables r^j for generic events simulated by PYTHIA 
and pure DCC events simulated with the 'snowball' model. Comparisions 
with the Ti/s obtained with binomially distributed pions and the l/(2y/f) 
classical limit of DCC's. 







PYTHIA 


DCC 


binomial 


V(2V7) 


% 


3 


r i,j i °Y M 


r i,j -'- °V»,7 






1 


1 


1.00 ± 0.02 


0.57 ± 0.01 


1.00 


0.50 


2 


1 


1.00 ± 0.05 


0.43 ± 0.03 


1.00 


0.33 


3 


1 


1.04 ± 0.13 


0.38 ± 0.05 


1.00 


0.25 





2 


1.36 ± 0.04 


1.55 ± 0.06 


1.36 


1.80 


1 


2 


1.36 ± 0.10 


0.66 ± 0.06 


1.30 


0.62 


2 


2 


1.47 ± 0.26 


0.44 ± 0.09 


1.25 


0.31 





3 


2.13 ± 0.25 


2.98 ± 0.39 


1.89 


3.54 


1 


3 


2.03 ± 0.43 


1.19 ± 0.31 


1.74 


0.90 





4 


3.06 ± 0.94 


6.82 ± 2.18 


2.70 


7.34 
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Table 2. The effect on the r^i of an admixture of DCC and generic 
(PYTHIA) events. DCC domains from the DCC-generator/GEANT are 
added to various fractions of random PYTHIA/GEANT events. The first 
column represents the fraction of events in which a DCC is overlaying a 



DCC fraction 


7*1,1 ± cr ri l 


r 2 ,i ± ov 21 


r 3; i ± a r31 


events 


0.00 


1.01 ± 0.02 


1.02 ± 0.05 


1.09 ± 0.14 


51741 


0.02 


1.00 ± 0.02 


1.00 ± 0.05 


1.01 ± 0.15 


51741 


0.05 


0.97 ± 0.02 


0.93 ± 0.05 


0.95 ± 0.10 


51741 


0.10 


0.95 ± 0.02 


0.89 ± 0.04 


0.89 ± 0.08 


51741 


0.20 


0.93 ± 0.02 


0.83 ± 0.04 


0.77 ± 0.07 


51741 


0.50 


0.84 ± 0.01 


0.71 ± 0.03 


0.68 ± 0.06 


40000 


1.00 


0.74 ± 0.01 


0.60 ± 0.03 


0.55 ± 0.06 


20000 
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